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Abstract 

The J-PARC E56 experiment aims to search for sterile neutrinos at the 
J-PARC Materials and Life Science Experimental Facility (MLE). In order to 
examine the feasibility of the experiment, we measured the background rates 
of different detector candidate sites, which are located at the third floor of 
the MLE, using a detector consisting of plastic scintillators with a fiducial 
mass of 500 kg. The gammas and neutrons induced by the beam as well as the 
backgrounds from the cosmic rays were measured, and the results are described 
in this article. 


Now at Quintessa K.K., Kanagawa, JAPAN 



1 Introduction 


We proposed a definite search for the existence of nentrino oscillations 
with Am^ near 1 eV^ at the J-PARC Materials and Life Science Exper¬ 
imental Facility (MLF) [1]. With the 3 GeV Rapid Cycling Synchrotron 
(RCS) and spallation neutron target, an intense neutrino beam from muon 
decay at rest (fiDAR) is available. Neutrinos come predominantly from 
the /i’*' decay : /ie"*" -|- -|- z/g. The oscillation to be searched for 

is —)■ Pe which is detected by the inverse /3 decay (IBD) interaction 

Pe -|- p —>■ -f n, followed by gammas from neutron capture. Figure 

(left) shows a bird eyes’ view of the MLF building. 
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Figure 1: A bird eyes’ view of the the MLF facility in J-PARC (left), and the 
measurement locations (red circles) of the experiment, (right; those are written as 
“Point 1”, “Point 2” and “Point 3”)) 

The unique features of the proposed experiment, compared with the 
prior LSND experiment [2] and experiments using conventional horn fo¬ 
cused beams (e.g.: m , are; 

1. The pulsed proton beam with about 600 ns spill width from J-PARC 
RCS and muon long lifetime allows us to select neutrinos from fxDAR 
(The protons are produced with a repetition rate of 25 Hz, where each 
spill contains two 100 ns wide pulses of protons spaced 600 ns apart). 

This can be easily achieved by gating out for about 1 ps from the 
start of the proton beam spill, that eliminates neutrinos from pion 
and kaon decay-in-fiight. 
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2. The time gate width with ~10 /zs to obtain the neutrinos from ^DAR 
can also provide superb rejection capability on the cosmic ray back¬ 
ground by a factor of ~ 4000. 

3. The spallation neutron source is the mercury target, which is high-Z 
material, surrounded by thick iron and concrete shields as shown in 
Fig-i Due to a strong nuclear absorption of vr and /r in the mer¬ 
cury target, neutrinos from decay are strongly suppressed up to 
about the 10“^ level. The resulting neutrino beam is predominantly 
Ue and from jj,'^ with contamination from other neutrino species 
at the level of 10“^. 



Target vessel M 
(Triple-walled structure) 


(a) The mercury target 


(b) Structure around the mercury target 


Figure 2: Mercury target (left) and the structure, which surrounds the target (right). 

4. Ue interacts via IBD and its cross section is known to a few percent 
accuracy |1]. 

5. The neutrino energy can be reconstructed from the positron visible 
energy by adding ~ 0.8 MeV. 

6. The Ufj_ and z/g fluxes have different and well dehned spectra. This 
allows us to separate z/g due to z/^ —)• Ue oscillations from those due 
to /i“ decay contamination. 

In addition to these, Gd-loaded liquid scintillator [5] is used to reduce 
the accidental backgrounds in the E56 experiment. Figure shows the 
current detector design. 

In the proposal, two detectors are placed on the third floor of MLF 
with a baseline ~20 meters from the mercury target. 

In order to examine the feasibility of the J-PARC E56 experiment, we 
carried out an on-site test experiment (MLF 2014BU1301 experiment), 
which was mainly dedicated to measure the beam related backgrounds. 

The required accuracy of the measurements is a few ten % to check the 
feasibility of the E56. The data was taken from April to July 2014 using a 
500 kg plastic scintillator. The setup and calibration of the detector, and 
the results of background rate measurements are described in this article. 
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Figure 3: Current design of the E56 25 tons detector. Left: schematic view; right 
3D drawings. 


2 Definition of the IBD Signal 

Before describing the results, the dehnition of the IBD signal inside the 
liquid scintillator is clarihed since it will be used often within the text. 
As explained before, the IBD interaction is z/g + p —)■ e"*" + n, followed 
by gammas from the neutron capture. For the detection of the positron, 
which we call “prompt signal region”, the energy is selected to be 20< 
E <60 MeV and the hit time of the activity inside the detector to be 
1< T <10 /is, where T is the time from the proton beam start timing. 
These criteria are based on the features of the fiDAR neutrinos. On the 
other hand, to catch the delayed gamma (called “delayed signal region”), 
we select the activity, which have the energy in the range 7< E <12 MeV 
and the hit time to be < 100 ps. 

To examine the background rate, these selection criteria are used unless 
noted otherwise. 


3 Plastic Scintillators with 500 kg Fiducial 
Volume (the 500kg Detector) 

The 500 kg plastic scintillator detector consists of the main target 
scintillators and two layers of charged particle vetoing system. We describe 
the setup, the calibration and the obtained resolution of the scintillator 
detector in the following subsections. 

3.1 Setup 

Figure l^shows a schematic view of the 500 kg plastic scintillator coun¬ 
ters placed at the third floor of MLF. The yellow parts of the Fig. |^show 
the 500 kg main scintillators in this experiment. They consist of two types 
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of scintillators: 12 pieces of 11.7/13.7 (trapezoid) x7.6 x 182 cm^ scintilla¬ 
tor (ID) and 12 pieces of 16.9/18.8 (trapezoid) x7.6x 182 cm^ scintillator 
(3D). To minimize the dead space, 3D scintillators were located on the 
both sides, while ID scintillators were located in the central part. Each 
end of the scintillators was viewed by two PMTs. Signals from each PMT 
were recorded by FINESSE 500 MHz 8 bit FADC [B]. The 500 kg detec¬ 
tor was surrounded by two layers of charged particle vetoing systems, the 
Inner and Outer vetoes. The Inner veto (red parts in Fig. covers the 
surfaces of the top, bottom and both sides of the 500 kg detector. The 
thickness of the plastic scintillator for Inner veto is 4.3 cm. The Outer 
veto (green parts in Fig. surrounds the 500 kg detector and Inner veto, 
with mostly 2 layers of 6-8 mm thick plastic scintillators. PMT signals 
from the veto counters were recorded by FINESSE 65 MHz FADC with 
50 ns RC-£Iter |6]. 



Figure 4: Schematic view of the 500 kg plastic scintillator detector at the third floor 
of MLF(left: front view, right: side view). The 500 kg plastic scintillators (yellow) 
were surrounded by two layers of charged particle vetoing system. Inner veto (red) 
and Outer veto(green). 


3.2 Calibration 

We used cosmic muons to calibrate the energy and timing and to mea¬ 
sure the attenuation length of the scintillator. Four pairs of plastic scin¬ 
tillators were prepared to dehne the cosmic ray tracks. Figure shows 
a schematic view of the cosmic muon trigger counters. The size of the 
scintillators located at both sides are 132 cm (length) x 10.5 cm (width) 
X 4.0 cm (thickness), while those located in the center are typically 80 cm 
(length) X 5.0 cm (width) x 1.5 cm (thickness). 

To measure the attenuation length of the scintillators, we made some 
dedicated runs in which we changed the position of the trigger counters. 
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Figure 5: Schematic view of the cosmic muon trigger counters. Four pairs of plastic 
scintillators were prepared to trigger cosmic muons. 


Figure shows the typical attenuation curves measured for each scin¬ 
tillator type. We measured the attenuation curve and parametrized for 
each scintillator. By considering the attenuation length, the reconstructed 
charge becomes independent from the incident position. 



Figure 6: Typical attenuation curve for each scintillation type (ID and 3D). We 
measured the curves and parametrized for each scintillator (magenta line. The £t 
function is the 2nd polynomial -|- exponential). 


The gain of each PMT was also calibrated by using cosmic ray events. 
The output from each scintillator is stable within 1-2% during the mea¬ 
surements with the gain correction. 

We also adjusted the timing of each PMT. Time offsets were deter¬ 
mined to minimize the time difference between each pair of PMTs on 
each end. The velocity of cosmic muons passing through the detector was 
considered, where the typical light velocity inside the scintillators is 14.3 
cm/ns. 
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3.3 Resolutions 

We also evaluated the position and energy resolutions of the 500 kg 
detector using the cosmic ray events. As the results of the precise timing 
calibrations, the obtained position resolution is = 2.6 cm for MIP 
energy. The typical energy resolutions of the scintillators are 3.3% for ID 
and 4.5% for 3D at the middle of the scintillators for MIP energy, which 
is about 13 MeV. 

3.4 Detector simulation 

The resolutions described above and other detector responses, such as 
quenching effects with Birks’ law, light attenuation and the threshold 
effect, were implemented to Geant4 [7] based Monte Carlo simulation 
(MC). 


3.5 Veto Efficiency 

In order to measure the efficiency of the Inner Veto (IV) and Outer 
Veto (OV) systems on the background measurement, their particle tagging 
efficiency was measured. The veto efficiency (e) is dehned as follows: 


N of coincident event OV, IV and Target (Triple Coincidence) 
N of coincident event OV and Target (Double Coincidence) 

N of coincident event OV, IV and Target (Triple Coincidence) 
N of coincident event IV and Target (Double Coincidence) 

The energy spectra of the main 500 kg detector for different veto con¬ 
ditions in no beam period are displayed in Fig Around the 80 MeV 
region of the energy distribution without veto, there is a peak created by 
the cosmic muons. Finally, the veto efficiency for different energy ranges is 
summarized in Table [T| where it can be seen that a total veto efficiency is 
better than 99.8% assuming the IV and OV inefficiencies are independent. 
In these Table and Figure, it is also possible to see that muons deposit an 
energy larger than 60 MeV and most of cosmic ray muon background are 
thus rejected. 


Table 1: Summary of IV and OV efficiency for different energy ranges. 


Energy Range [MeV] 

Ov 

^ov 

(IV or OV) 

20 < E < 60 

96.8±0.2% 

94.1±0.2% 

-99.8% 

60 < E < 100 

99.5±0.04% 

96.2±0.1% 

—99.9% 

100 < E <140 

99.6±0.07% 

95.1±0.3% 

—99.9% 
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Figure 7: Energy spectra for the 500 kg detector for different veto conditions. The 
peak due to the cosmic muons are seen around 80 MeV for the distribution without 
veto. 


4 Background Measurement Locations 

Figure shows the overview of the MLF building (left), and the mea¬ 
surement locations (right). The 500 kg detector was moved to measure 
the backgrounds at each point. The baselines from the mercury target are 
~17 m, ~20 m, ~34 m for “Point 1”, “2” and “3” in Fig. respectively. 
We accumulated data for two weeks per each point. 

Only the results of “Point 2” and “Point 3” are described in this article. 
The results of “Point 1” is expected to be published later since the analysis 
is complicated due to large amount of the neutron background. 

The number of spills during the beam to be used for the analysis are 
19,545,739 for “Point 2” and 26,225,816 for “Point 3”, respectively. 

5 Fast Neutrons Background from Beam 

One of the main purposes of the background measurement is to mea¬ 
sure the Michel electron background induced by beam fast neutrons, which 
was indicated by the previous background measurement at the hrst floor 
of MLF building [1]. According to Geant4, fast neutrons whose kinetic 
energy are larger than 200 MeV can produce charged pions. The charged 
pions create the Michel positrons, which has the same energy and hit tim¬ 
ing as IBD events in the vr —)■ p —>• e decay-chain and the scattered fast 
neutrons are thermalized and captured by Gd. Therefore, this could be a 
serious correlated background for the E56 experiment. The flux of such 
fast neutrons at the third floor of MLF was estimated to be some orders 
of magnitude smaller than that in the hrst hoor [T]. We briehy describe 
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backgrounds signal 


the basic idea and backgrounds of this measurement. 

Figure]^ shows the dehnitions of the “signal” and “background” of this 
measurement. The signal is a Michel electron induced by beam fast neu¬ 
trons. Fast neutrons coming on the (proton) beam bunch timing hit our 
detector and produce pions. These pions then produce Michel electrons 
in the decay chain of n -|- p (or C) —)■ X + tt^, then /i+ —)■ e+. 

The signature of the signal is thus the coincidence between an activity 
without any veto hits on the bunch timing and a “prompt signal” about 
2.2 fis later from the beam timing. Backgrounds for this measurement are 
clipping cosmic muons, Michel electrons from cosmic muons and neutral 
particles (gammas and neutrons) from cosmic rays as shown in Fig. 
and they come either during beam-on or beam-off. On the other hand, 
signal comes only when the beam is on. The basic idea of this measure¬ 
ment is thus to extract signals from backgrounds by subtracting beam-off 
activities from beam-on activities. 
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Figure 8: The dehnition of “signal” and “backgrounds” of this measurement. Signal 
is a Michel electron induced by beam fast neutrons. Backgrounds are clipping cosmic 
muons, Michel electrons from cosmic muons and neutral particles from cosmic rays. 


Based on the concept of the measurement, we took the following data 

• beam-on: To observe activities on the beam bunch and around the 
bunch timing (the timing dehnition is described later); 

• beam-off (off-timing) : To subtract backgrounds from beam-on data 
We took data 20 ms after each beam bunch spill. Because detector 
responses such as the PMT gains, efficiency of veto counters and oth¬ 
ers are exactly the same with the last beam spill, we can subtract the 
backgrounds from the beam signals without systematic uncertainties. 

• no-beam: When the accelerator is off. To evaluate truly beam unre¬ 
lated background; 
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• cosmic muons: To calibrate the detector. 

Note that the radio frequency (RF) of the RCS provides the beam timing 
signal for the measurement, thus we used it to take data around the beam 
and 20 ms later. For the no-beam data, the clock trigger was used. 


5.1 Measurement 


The search for beam neutrons which induce Michel electrons was per¬ 
formed by detecting their prompt signals. Figure]^ (a) and (b) show the 
correlation between the energy and the timing of the events observed at 
“Point 2” and “Point 3”. The zero of the horizontal axis corresponds 
to the starting time of the proton beam. We observed the large number 
of activities around the beam time, however the number of activities are 
decreased rapidly as a function of time. Also no-beam data is shown in 
Fig. I (c) for comparison. Here we used the clock trigger, thus the time 
axis is arbitrary. In this measurement the prompt signal is dehned by: 

• 20 < E[MeV] < 60 

• 1.75 < f[/is] < 4.65 from the rising edge of the hrst beam bunch. 


to estimate the number of activities in the IBD prompt region. 

As described before, we compare the beam-on and beam-off data to 
subtract other activities without any systematic uncertainties. A huge 
number of the clipping muon background (order of a few 100 Hz) was 
rejected by applying the charged veto cut. Figures IT show the energy 
distributions of events in the prompt timing window, 1.75 < f[/rs] < 4.65 
from the beam bunches, and beam-off data, before and after applying 
the charged veto cut. The observed rates are summarized in Table 
The numbers of events in the prompt energy range are both consistent 
statistically between beam-on and beam-off data either with or without 
applying the charged veto cut. 


Table 2: Summary of the background rates in the prompt region. 



“Point 2” (xlO-Vspill) 

“Point 3” (xlO-Vspill) 

beam-on w/o veto 

16.8 ±0.3 

15.3 ±0.2 

beam-off w/o veto 

16.4 ±0.3 

15.4 ±0.2 

subtraction (on-off) 

0.4±0.4 

-0.1±0.3 

beam-on w/ veto 

1.58 ±0.09 

1.41 ±0.07 

beam-off w veto 

1.52 ±0.09 

1.33 ±0.07 

subtraction (on-off) 

0.06±0.13 

0.08±0.10 


To improve the sensitivity, an additional cut was applied before ob¬ 
taining the final result. Figure IT shows a simulated energy distribution 
on the bunch timing when fast neutrons produce charged pions in the 500 
kg detector. We can expect some events on the bunch timing associated 
with beam Michel electron backgrounds. On the other hand, as shown in 
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(c) “no-beam” 

Figure 9: Correlation between energy and timing of the events observed at 
(a) “Point 2”, (b) “Point 3” and (c) “no-beam”. Note that the vertical axis of the 
two dimensional plot is based on the log-scale, and one bin of the vertical axis is 
divided equally in the log-scale. 
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(b) “Point 3” 

Figure 10: Energy distributions of events in the prompt timing window and beam-off 
data, before and after applying the charged veto cut. (a) “Point 2” and (b) “Point 3”. 
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Figure 11: Energy deposit on the bunch timing for beam Michel electron with the 
neutron kinetic energy 300 to 500 MeV (flat). Note this is MC. Red, blue and black 
lines are the energy deposit from neutron interactions, pion decays and their sum in 
each event. Note that one bin of the horizontal axis is divided equally in the log- 
scale, and the unit of the vertical axis is arbitrary. A peak around 3 MeV in the 
blue histogram corresponds to the muon kinetic energy, 4.2 MeV, from stopped pion 
decay including the quenching effects with the Birks’ law. The quenching effects, 
light attenuation in the scintillator and other detector responses such as resolutions 


and threshold effects are implemented to the Geant4 as described in Section 3.4 
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Fig. the beam unrelated backgrounds have no activities on the bunch 
timing. Michel electrons from cosmic muons can have an activity on the 
bunch timing only when the parent muon comes on the bunch timing ac¬ 
cidentally. However, since the muon is a charged particle, it can be easily 
rejected by the veto counters. We can thus strongly suppress backgrounds 
by requiring on-bunch activities without hits in the veto counters, and at 
least one on-bunch activity (Edep>4 MeV) with this veto conditions was 


required. Figure 1^ shows the estimated selection efficiency of this on- 
bunch cut as a function of the incident neutron kinetic energy based on 
MC. Though most of the events have more than 4 MeV energy deposit 
at on the bunch timing, a part of the events are rejected by self-vetoing. 
The selection efficiency has slight dependence on the incident neutron ki¬ 
netic energy. Because we do not know the energy spectrum of the incident 
neutron well, we assumed the selection efficiency, Conbunch = 0.9. 
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Figure 12: Estimated selection efficiency of the on-bunch cut as a function of the 
incident neutron kinetic energy. 


Figure [T^ shows the energy distributions after applying the cut to the 
500 kg detector data. The beam unrelated distribution were obtained 
with the energy distribution of beam-off data without applying the on- 
bunch cut, and the accidental coincidence probability: hit rate of neutral 
activities on the bunch timing. The mean hit rates without veto activ¬ 
ities on the bunch timing were 3.1% (“Point 2”) and 0.6% (“Point 3”). 
The observed event rates during beam-on were (4.60 ± 1.53) x 10“^/spill 
(“Point 2”) and (1.53±0.76) x 10“^/spill (“Point 3”), while the estimated 
event rate by beam unrelated activities were (4.91 ± 0.28) x 10“'^/spill 
(“Point 2”) and (0.86 ± 0.04) x 10“’^/spill (“Point 3”), and both rates 
between the observed and predicted numbers are consistent. By consid¬ 
ering the efficiency of the on-bunch cut, eonbunch = 0.9, the upper limit of 
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the event rate of the beam Michel electron are thus 2.5 x 10“^/spill (90% 
C.L., “Point 2”) and 2.1 x 10“^/spill (“Point 3”). Table summarizes 
these numbers. 



energy[MeV] 


(a) “Point 2” 



energy[MeV] 


(b) “Point 3” 

Figure 13: Energy distributions after applying the on-bunch cut. The beam unre¬ 
lated distribution were obtained with the accidental coincidence probability, 3.1% 
(“Point 2”) and 0.6% (“Point 3”), and the energy distribution of beam-off data with¬ 
out applying the on-bunch cut. 


6 Accidental Background 

Another important background comes from accidental coincidences. 
The accidental background rate is estimated using the multiplication of 
the single rates of the background in the prompt region and that in the de- 
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Table 3: Summary of the background rates in the prompt region after the on-bunch 
cut._ 



“Point 2” (xlO-Vspill) 

“Point 3” (xlO-Vspill) 

Beam-on data 

4.60 ± 1.53 

1.53 ±0.76 

Prediction 

4.91 ±0.28 

0.86 ±0.04 

Subtraction 

-0.31T1.56 

0.67±0.76 

90% C.L. upper limit 

< 2.5 

< 2.1 


layed region. Therefore, an absolute rate of the background measurements 
on the prompt and the delayed region with the 500 kg detector are crucial. 
Measurements with small size detectors are also important to estimate the 
contents (particle identihcation: PID) of the prompt background. 

We hrst explain the background measurements for the prompt region, 
and then for the delayed region. 


6.1 Single Background Rate Measurement for Prompt 
Region 

As shown in the previous section, the background rates in the prompt 
region after charged particle veto are consistent statistically between “beam- 
on” and “beam-off” data. This means that the beam related background 
is small enough compared to the beam unrelated background which is 
expected to be mainly due to the neutral particles, gammas or neutrons 
induced by the cosmic rays. Since the response of the E56 detector is quite 
different for gammas and neutrons, it is essential to estimate their influ¬ 
ences separately in the prompt region. Thus measurements of the particle 
identihcation for the prompt region was done at Tohoku University using 
small size detectors at hrst, which were used to construct and to check 
the hux model. Finally, the hux model was compared with the 500 kg 
detector data taken at the third hoor of MLF. All details are described in 
the following subsections. 


6.1.1 PID and Energy Measurements with Small Size Detectors 


The calibrations of the small size detectors are described elsewhere [8] , 
so only a summary of the measurements is shown here. 

There are no reliable gamma hux measurements or models so far at 
the ground level, which contributes to the prompt event of the acciden¬ 
tal background. The high energy gammas are supposed to be produced 
mainly by the decay of which are created by /x-nucleon interactions in 
the walls of the building or air, outside the detectors. Therefore we con¬ 
structed an empirical gamma hux model with very simple two exponential 
functions for energy at the surface of the detector at hrst. 
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Nal Energy 



Figure 14: Nal energy spectrum for the remaining events after applying the veto cut. 
The environmental gammas below 3 MeV are generated based on the citation [U]. Note 
that the horizontal axis is based on the log-scale, but one bin of the axis is divided 
by 0.5 MeV. 


Parameters, A, B, C, D were determined from a Sodium Iodide (Nal) (a 
cylinder with 2” diameter and 2” height) detector data at Tohoku. Cos¬ 
mic ray muons were rejected using plastic scintillator anti-counters, which 
surround the Nal counters, with a veto efficiency better than 99%. Note 
that the cosmic fi veto also eliminated the gammas from tt^s created by /i- 
nucleon interactions inside the detector. This statement can be expanded 
to any size detectors in general. Figure 14 shows the results of the mea¬ 
surement, where the energy spectrum in the yellow part can be described 
by two exponential functions with mean energies of 3 (“A” in Eq. and 
26 (“5”) MeV (the MC gammas are generated using the two exponential 
functions from the Nal surface assuming an isotropic direction). Com¬ 
paring data and MC spectra, these components rates are measured to be 
150 (“C”) for the first exponential function and 25 (“H”) Hz/m^ for the 
second, respectively. 

Another important measurement had been done using a small liquid 
scintillator (NE213 110]) which can separate neutrons from gammas effi¬ 
ciently using a pulse shape of the scintillation light (PSD; Pulse Shape 
Discrimination). Neutron signals have larger tails than those from gam¬ 
mas, thus we can check the flux of fast neutron in addition to the gamma 
flux which was measured by Nal already. Fast neutrons induced by cosmic 
rays are measured and modeled by some previous works. In this work we 
employed the model and the parameters given by Wang et al [T2|, who 
gave the empirical functions on the kinetic energy (E„), the multiplicity 
(M) and the zenith angle {9) of the fast neutrons as a function of the 
muon energy {E^). Figure 15 shows those distributions. 
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Figure 15: Distributions of the neutron kinetic energy (upper left plot), multiplicity 
(upper right plot) and zenith angle (lower plot) of the cosmic muon induced fast 
neutron events generated by the empirical functions in the reference [12]. 


A cylindrical aluminum housing (5” of diameter and 2” height), with 
white painted inner walls, is hlled with NE213, closed with a glass plate 
and attached to a 5” PMT (R1250-03). The NE213 detector was also 
surrounded by the plastic scintillators, and Fig. [T^ shows the results of 
the measurement. 

With this set-up, the fast neutron flux detected above 20 MeV was of 
(1.28 ± 0.05) X 10“^ Hz (statistical uncertainty), while the MC based on 
the reference na gives a rate of 1.12 x 10 ^ Hz. For gammas, in the same 
energy range, the measured rate was of (1.18 ± 0.04) x 10“^ Hz (statisti¬ 
cal uncertainty), while the MC based on the Nal measurement described 
above gives 0.95 x 10“^ Hz. MC simulation generates the gammas at 
the surfaces of the detectors isotropically. Therefore, the data and MC 
above 20 MeV for both gammas and neutrons agrees within 20% of the 
uncertainty. 


6.1.2 Measurement for Prompt Background with 500 kg De¬ 
tector 

Figure [T^ shows the energy distribution of no-beam data after applying 
charged particle veto and the Michel electron cuts. Because the activities 
induced by cosmic rays are dominant in the prompt region as described 
above, no beam data is used to estimate the background for the simplicity. 
The MC estimation with cuts is overlaid. Fast neutrons are generated 
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Figure 16: Reconstructed energy distribution for events with low or no energy 
deposition in the veto. The data and MC components of neutrons and gammas are 
also compared. Data events without veto applied are also shown for comparison. 


based on the model in the reference [12], and the gammas are generated 
based on the Eq.Q at the surface of the 500 kg detector isotropically 
again. 

The cosmic muons were eliminated by applying the charged particle 
veto cut because of the high veto efficiency (>99.8%) as described in Sec¬ 
tion |3.5[ The Michel electrons can be rejected by detecting the parent 


muons coming earlier than the prompt events by 70 ns. Figure 17 shows 
the timing distribution of parent muons, which generate Michel electrons 
within a time window, 3.9 < t[fis\ < 5.4, based on a toy MC. This time 
window is applied to the real data observation to obtain Fig. [l^ Some 


of the events within the time window remained after the cut because the 
Michel electrons came too close to their parent muons (corresponding to 
the right red hatched histogram of Fig. 17 ) Events in the negative 

region also remained because they were out of FADC time window (cor¬ 
responding to the left red hatched histogram of Fig. mi' The rejection 
power of the Michel electron cut was thus estimated to be 7. After the 
cuts, the remaining events are composed of cosmic gammas and neutrons 
mainly for the prompt region. 

The cosmic induced fast neutrons make the correlated background, 
however the 25 ton detectors for the E56 will have PID capability, which 
can reject fast neutrons relative to electrons and gammas by a factor of 
more than 100 [1]. Neutrons in the remained events are thus not harmful 
and only gammas can be a prompt background. By using gamma flux 
model described above, the number of gammas in the remained events 
was equivalent to (6.4 ± 0.5) x 10“® /2.5/is. 

Note that the predicted rate from the measurement at Tohoku Uni- 
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Figure 17: Timing distribution of parent muons, which generate Michel electrons 
within a time window, obtained by a toy Monte Carlo. The time window for Michel 
electrons is 3.9 < t[/rs] < 5.4. 
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Figure 18: Energy distribution of no-beam data after applying charged particle veto 
cut and the cosmic Michel electron cut. 
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Figure 19: 


Self shield effect of liquid scintillator as function of the shield thickness. 


versity was consistent with the rate measured by the 500 kg detector at 
the MLF within 6%. As mentioned previously, the goal of the test exper¬ 
iment is to check the energy and the rate of the on-site backgrounds. The 
neutron and gamma flux models in this work can be used to discuss the 
feasibility of the E56 experiment |8]. 

6.2 Single Background Rate Measurement for De¬ 
layed Region 

There are two dominant backgrounds for the delayed region. They are 
neutrons and gammas induced by the beam, and are described in the next 
subsections. 

6.2.1 Beam Neutrons 

Thermal neutrons in the third floor of MLF is not harmful for the E56 
experiment since they are stopped at the buffer region of the E56 detector, 
whose thickness of liquid scintillator is 50 cm (Fig. |^. But the middle 
energy (10 MeV < Ekin) neutrons on the beam bunch can be thermalized 
inside the detector, and can be captured in the liquid scintillator and emits 
the gammas inside the hducial volume. They become backgrounds for the 
delayed region. Figure shows the attenuation rates of thermal neutrons 
inside the liquid scintillator generated by Geant4. Most of neutrons which 
have less than 10 MeV are rejected at the buffer region. 

We estimated the neutron flux on the beam bunch timing above 15 
MeV using data assuming the reaction rates of neutrons inside the 500 kg 
detector based on Geant4 since the beam gammas are dominated below 15 
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MeV as described later. Here we used the number of activities and their 
energy spectrum with the time window less than 3 /rs from the beam start 
time to estimate the flux. The MC simulation also assumed that such fast 
neutrons are created at the mercury target isotropically and comes to 
the 500 kg detector directly. Figure 20 shows the schematic view of the 
assumption. 



Figure 20: Geometry of the mercury target and the 500 kg detector. 


This assumption is valid for the neutron backgrounds with only similar 
solid angle positions since the scattered angle, energy and production rate 
are tightly correlated, however the assumption has a good shape for the 
“Point 2” and “Point 3” as described later. 

The flux, 0(i?n), is denoted as 


0(^n) = ^exp(-En/30), 


( 2 ) 


where is the neutron kinetic energy and a is a number of neutrons 
that are generated at the mercury target in a spill. Note that a is just an 
“effective” number of generated neutrons at the target since the 500 kg 
detector is put on the third floor of MLF, and the neutrons are highly 
attenuated by the iron or concrete, which surrounds the target as shown 
in Fig. a is estimated to be (387±12)/spill as a £t result at “Point 2”. 
As shown in Fig. this flux reproduces well the measurement above 15 
MeV, which is dominated by activities from neutrons. 

The neutron rate at “Point 2” with the energy above 15 MeV is 


(5.4±0.02) X10 ^/spill/0.3MW/500kg, while that at “Point 3” is 
(2.0±0.01) X 10“^/spill/0.3MW/500kg. Although the neutron flux depends 
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Figure 21: Energy distribution of the 500 kg detector for on-bunch timing. The 
black cross is the measured energy spectrum and the red one is a MC htted to the 
measured one. 


on the scattering angle from the mercury target, these are well explained 
by the 1/r^ law, where r is the detector baseline, within 5%. 

6.2.2 Beam Gammas 

Beam gammas are produced by capture reactions of the thermalized 
beam neutrons at the materials of MLF building such as concrete walls or 
floors. They are injected into the E56 detector, therefore they are crucial 
background in the delayed region. 

As shown in Fig. [^a) and (b), the hit time distributions of the low 
energy activities were almost flat. Therefore even the rate measurements 
during short time span such as only 1 /rs is useful to estimate the back¬ 
ground rates for the IBD delayed region. Figure shows the energy spec¬ 
tra below 16 MeV measured by the 500 kg detector. 4.35 < t[fis] < 5.35 
from the rising edge of the hrst beam bunch shown in Fig. [^a) and 
(b) were used to estimate background in the IBD delayed region. In 
the energy region of 7 < i7[MeV] < 12, the background rate with the 
500 kg detector are (8.26 ± 0.07) x 10“^//is/0.3MW at “Point 2”, and 
(3.65 ± 0.12) X 10-V/rs/0.3MW at “Point 3”. 

As shown in the bottom plot of Fig. [^ there is a room for the re¬ 
mote maintenance of the mercury target under “Point 2” and “Point 3”. 
The floor of “Point 2” is made of concrete which has 120 cm thickness 
and is thinner than that of “Point 3” by 30 cm. This structure of the 
MLF building can make the difference of the background rates between 
“Point 2” and “Point 3”. 
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Figure 22: Energy distribution of the backgrounds for delayed region measured with 
the 500 kg detector. Black: energy spectrum at “Point 2”, red: that at “Point 3”, 
blue: no beam data. 



Figure 23: The top view of the third floor of MLF around the measurement points. 
Bottom is the side view. The floor of “Point 2” is made of concrete which has 120 
cm thickness and is thinner than that of “Point 3” by 30 cm. 
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7 Summary 

The J-PARC MLF 2013BU1013 test experiment was carried out to 
measure the backgrounds and to examine the feasibility of the J-PARC 
E56 experiment. The results show that: 

• Beam fast neutrons: 

The beam fast neutrons which create charged pions inside a 500 kg 
plastic scintillators, are not observed during 2 weeks. The 90% C.L. 
upper limit is set to be 2.5xl0“^/2.9/is/0.3MW at “Point 2”, and 
2.1xl0-V2.9/is/0.3MW at “Point 3”. 

• Accidental: 

1. The gammas or neutrons induced by cosmic rays are observed for 
the prompt region. The neutron rate is the same as expected in 
the reference [12] . Gammas are newly recognized by this experi¬ 
ment. The number of gammas is equivalent to (6.4 ±0.5) x 10“® 
/2.5/rs in the 500 kg detector in the prompt region. 

2. Neutrons and gammas induced by the beam are observed for 
the delayed region by the 500 kg detector. The flux of neutrons 
are estimated by the data, and the shape of the flux is consis¬ 
tent with the exponential with a slope constant of 30 MeV. The 
effective neutron production rate is (387±12)/spill at the tar¬ 
get. “Point 2” and “Point 3” have same neutron flux within 
5%, considering the baseline effect. The rates of gammas are 
(8.26 ±0.07) X 10-V/is/0-3MW (“Point 2”) and (3.65 ±0.12) x 
10“®//is/0.3MW (“Point 3”) in the energy region of 7 < i7[MeV] < 
12 . 

These information have been translated to the background rates at the 
E56 detectors using the MC simulation with evaluation of the self-shielding 
effects, and the results show that the backgrounds at the “Point 2” provide 
no specihc problems to perform the J-PARC E56 experiment [S]. 
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